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SUMMARY 

Pulsed  Detonation  Engines  (PDFs)  cover  an  unusually  wide  range  of  physics  and  thermochemistry  . 
1  he  flow  is  time  dependent  and  compressible  and  the  chemical  structure  of  the  fuel  has  been  shown 
experimentally  to  have  a  leading  order  influence  on  the  transition  process  Topics  of  direct  relevance 
to  device  design  include  fuel  mixture  injection  and  obstacle  enhanced  localised  turbulent  explosions. 
The  latter  should  result  in  a  transition  to  detonation  and  propagation  of  the  resulting  wave.  Calculation 
methods  -  though  more  rarely  adequate  chemistry  descriptions  exist  for  laminar  shockwave 
containing  flow-s  and  the  work  performed  here  is  aimed  to  study  events  leading  up  to  the  crucial  DDT 
phase  During  the  latter  a  significant  traction  of  the  total  cycle  time  may  be  expended  with  the  less 
reactive  fael-3ir  mixtures  targeted  in  practical  applications.  Work  to  date  has  show  n  that  the  creation 
of  the  initial  explosion  kernel  is  strongly  dependent  upon  the  state  of  the  mixture  at  ignition  and  the 
duration  of  the  injection  phase.  Work  performed  in  an  earlier  contract  F6i  775-00- WE054  has  shown 
that  limiting  expressions  in  the  high  Damkohler  number  regime,  combined  with  closures  at  the  second 
moment  level  for  velocity  and  scalar  fields,  can  reproduce  comparatively  detailed  experimental  data 
w  ith  reasonable  accuracy.  It  was  also  noted  that  in  the  key  initial  phase  the  alternative  limiting  source 
term  forms  essentially  bracketed  experimentally  observed  trends.  Work  on  the  same  contract  also 
showed  that  large  scalar  space  transported  PDF  approaches  can  capture  key  aspects  of  premixed 
parabolic  turbulent  flames  at  high  Reynolds  numbers. 

The  w  ork  reported  here,  as  pan  of  the  current  contract  (F61775  -  01  WE050)  extends  the  past 
effort  in  two  directions.  The  first  is  related  to  the  sensitivity  of  the  initial  explosion  phase  (prior  to 
DDT)  to  the  state  of  the  mixture  resulting  from  injection  of  the  reactant  mixture.  The  work  naturally 
highlights  the  need  to  reconcile  the  .uniting  source  term  closures  controlling  the  heat  release  as 
outlined  in  previous  work.  The  second  aspect  of  the  current  work  features  computations  of  two- 
dimensional  elliptic  unsteady  flows  w  ith  comprehensive  chemistry  and  a  transported  PDF  approach 
closed  at  the  joint  scalar  level.  This  part  of  the  w'ork  is  exceptionally  resource  intensive  and 
constitutes  an  essential  "proof  of  concept"  study  aimed  at  providing  a  computational  procedure  in 
which  the  chemistry  of  the  fuel  appears  in  a  closed  form.  The  study  requires  significant  computational 
resources  and  the  test  case  has  therefore  been  chosen  to  draw  on  related  work  in  order  to  be 
computationally  tractable  within  the  current  framework  The  results  obtained  illustrate  that  the 
approach  is  technically  possible  and  constitutes  a  potential  route  for  dealing  with  the  chemical  kinetic 
effects  associated  w  ith  detonation  initiation  in  a  strongly  turbulent  flow  field.  However,  it  must  be 
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recognised  that  more  extensive  studies  are  required.  Permission  to  compute  the  actual  geometry  via 
ITAR  Export  Controlled  informaiion  permission  had  to  be  abandoned  in  favout  of  generic 
calculations  exploring  the  influence  of  turbulence  levels  on  the  strength  and  timing  of  the  onset  of  the 
DDT  process  The  results  show  a  strong  influence  and  indicate  a  possible  solution  to  the  problem  of 
detonation  initiation.  It  is.  however,  also  recognised  that  computations  of  the  sort  attempted  stretch 
our  current  understanding  of  turbulent  combustion  to  the  limit.  Accordingly,  experimental  studies  that 
can  further  highlight  the  physics  are  strongly  recommended. 

WORK  PROGRAMME 

TASK  1:  THE  SENSITIVITA  OF  THE  EXPLOSION  PHASE  TO  THE  STATE  OF 
THE  MIX  1 1  RE  RESULTING  FROM  INJECTION  OF  REACTANTS. 

Significant  new  work  aimed  at  reconciling  past  efforts  in  the  area  of  the  application  of  higher  moment 
based  methods  to  turbulent  gaseous  explosions  have  been  performed.  The  results  and  methodology 
are  outlined  in  Appendix  A  and  attached  to  the  current  final  report  The  enclosed  appendix  has 
appeared  in  CONFINED  DETONATIONS  AND  PULSED  DETONATION  ENGINES".  Editors 
CLD.  Roy.  S  M  Frolov.  RJ.  Santoro  and  S.A  Tsyganov.  Torus  Press.  Moscow.  2003.  The  final  task 
of  the  current  protect,  the  computation  of  geometries  of  direct  relevance  to  current  development  work 
on  PDF  engines,  had  to  be  deleted  due  to  an  absence  of  ITAR  Export  Controlled  information 
clearance.  Tnc  work  outlined  below  has  been  substituted  in  order  to  explore  the  influence  of  initial 
turbulence  on  the  strength  of  the  initial  DDT  kernel 

EFFECTS  OF  INITIAL  TURBULENCE  ON  THE  TIME  EVOLUTION  OF  EXPLOSION 
KERNELS 

'I  he  calculation  method  outlined  in  Appendix  A  features  closures  at  the  second  moment  level  for 
velocity  and  scalar  fields.  The  “return  to  isotropy"  and  “strain"  redistribution  parts  are  obtained  from 
model  formulations  derived  in  constant  density  flows.  I  he  model  used  here  is  that  of  Haworth  and 
Pope  (19S~i  which  is  based  on  a  stochastic  Lagrangian  Generalized  Langcvin  Model  (GLM).  The 
latter  features  a  non-linear  return  to  isotropy  and  has  been  calibrated  and  tested  for  a  range  o' 
homogeneous  turbulence  and  free  shear  flows.  An  important  feature  of  the  model  is  that  it  also  yields 
a  closure  for  the  pressure  scrambling  term  in  the  scalar  flux  equation  and  thus  allows  lor  realizable 
modelling  of  both  pressure  strain  correlations 
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However,  in  turbulent  premixed  flames,  the  "production"  by  strain  is  not  the  principal  cause  of 
anisotropy  and  it  has  been  pointed  out  by  Lindstedt  and  Vaos  (1999)  that,  regardless  of  specific  model 
formulations,  the  rate  at  which  energy'  is  transferred  between  components  is  substantially  lower  than 
the  rate  at  which  anisotropy  is  produced.  Instead,  the  dominant  influence  is  from  the  mean  pressure 
gradient  terms  <t>  and  associated  w  ith  preferential  acceleration  effects 


Deflagration  to  Detonation  Transition  Processes  in  Pulsed  Detonation  Engines 


The  potential  importance  of  terms  of  the  above  type  is  obvious  m  environments  with  strong  pressure 
gradients,  i  imhermore.  the  current  model  for  redistribution  should  also  include  the  contribution  from 
preferential  acceleration  effects.  The  functional  forms  arc  given  below 
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Ihe  isotropization  of  production  (IP)  model  is  adopted  for  the  closure  of  the  tensor  A.  The  model 
applied  here  :s  zeroth  order  in  the  anisotropy  tensor: 
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The  inclusion  of  terms  of  the  above  type  into  calculation  procedures  :>  r.n-tnviai  due  to  the  complex 
mathematical  forms.  The  fractal  based  model  of  Lindstedt  and  Vans  <1999)  has  been  extended  by 
Kuan  el  al  (2003),  as  shown  in  Appendix  I.  to  encompass  strongly  turbulent  flows  and  it  has  been 
shown  that  the  model  can  reproduce  explosion  kernels  with  overpressures  in  excess  of  160  kPa. 
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The  constants  C/  and  C:  have  been  calibrated  to  the  values  of  2.6  and  1.2  on  the  basis  of  flames 
stabilized  m  an  opposed  jet  geometry  (Lindstedt  and  Vans  1999).  The  velocity  reaction  rate 
correlation  term  is  closed  through  the  assumption  of  a  bimodal  PDF  and  the  thermochemistry  (eg. 
laminar  burning  velocities)  has  been  obtained  through  the  simulation  of  laminar  flames  using  detailed 
chemistry. 


RHSl  I  TS 

As  shown  in  Appendix  A.  good  agreement  can  be  obtained  with  experimental  data  for  the  base  case 
with  low  initial  turbulence  levels.  Increasing  the  initial  turbulence  levels,  achievable  in  practice 
through  a  variety  of  means,  has  a  strong  impact  upon  both  the  explosion  over-pressure  and  the  time 
taken  to  reach  this  point.  Ihe  predicted  pressure  traces  shown  in  Figures  1.1  and  1.2  corresponds  to 
initial  turbulence  fields  of  u”  1,  10  and  25  m  s. 

For  the  current  rather  weakly  reacting  stoichiometric  methane-air  mixtures,  the  time  to 
explosion  can  be  reduced  from  around  50  ms  to  10  ms  with  an  increase  in  the  turbulence  intensity 
from  1  m  s  to  25  m  s.  It  is  worth  observing  that  the  predominant  delay  is  caused  by  the  initial  flame 
travel  of  close  to  0.5  m  before  the  single  obstacle  is  reached  (see  Appendix  \  i  Turbulence  velocities 
of  the  order  2*  m  s  should  be  comparatively  readily  achieved  in  a  PDF  given  the  amount  of  flow  that 
can  be  generated  through  diverting  some  of  the  energy  of  each  detonation  pulse.  A  remaining 
principal  difficulty  is  to  assess  the  influence  of  turbulent  quenching  via  strain  and  direct  chemical 
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kinettc  effects  during  the  DDT  phase  under  such  conditions  and  further  experimental  and  supporting 
theoretical  work  remains  necessary. 


1  igure  1.1  Comparison  of  computed  and  measured  (dashed  line  l  pressure  traces  at  the  location  of  the 
obstacle  (0.4 — S  m)  for  turbulent  gaseous  explosions  in  pre-existing  turbulence  fields.  The  black  lines, 
red  lines  and  blue  lines  represent  calc  -  with  initial  lurbuk  3  is  of  u"  I.  10  and  25 
respectively 
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Figure  1.2  Comparison  of  computed  and  measured  (dashed  line)  pre-sure  traces  downstream  of  the 
obstacle  i  ^.208  m)  for  turbulent  gaseous  explosions  in  a  pre-existing  turbulence  fields.  Tbe  black 
lines,  red  lines  and  blue  lines  represent  calculations  with  initial  turbulence  fields  of  u”  =  1.  10  and  25 
m  s  respectively. 

The  corresponding  velocity  field  predictions  are  shown  in  Figures  1.3  to  1.6.  lhe  profile  shapes  and 
magnitudes  of  the  mean  velocities  illustrate  the  strong  influence  of  the  initial  turbulence  field  upon  the 
temporal  evolution.  The  very  strong  pressure  and  velocity  fluctuations  present  are  evident  both  in  the 
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pressure  traces  of  Figures  l  .l  and  1.2  and  in  the  mean  velocities  shown  in  Figures  1.3  and  14  The 
initial  dip  in  the  axial  velocity  profiles  correspond  to  the  peak  pressure  of  Figure  1.1  and  inaiks  the 
unset  of  explosion  With  higher  initial  turbulence  levels,  the  strength  of  the  explosion  increases  and 
this  is  reflected  in  the  magnitude  of  the  flow  acceleration  at  the  onset  of  the  explosion. 
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Figure  1.3  Comparison  of  computed  and  measured  (symbols)  axial  mean  velocities  at  the  location  of 
the  obstacle  for  turbulent  gaseous  explosions  in  pre-existing  turbulence  tields.  l  ines  are  as  in  Figure 
1.1 
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Figure  1.4  Comparison  of  computed  and  measured  (symbols)  cross-stream  mean  velocities  at  the 
location  of  the  obstacle  for  turbulent  gaseous  explosions  in  pre-existing  turbulence  fields.  Lines  are  as 
in  Figure  1.1. 

The  significant  increase  in  turbulence  levels  is  also  well  illustrated  in  Figures  1.5  and  1.6.  llowexer. 
significant  caution  is  required  in  interpreting  such  traces  given  the  extreme  complexity  of  the  How 
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field  Nevertheless,  it  is  worth  observing  that  for  the  case  where  experimental  data  is  available, 
surprisingly  good  agreement  is  obtained  as  shown  in  Appendix  A. 
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figure  1.5  Comparison  of  computed  and  measured  (symbols)  axis:  turbulence  velocities  at  the 
location  of  the  obstacle  for  turbulent  gaseous  explosions  in  pre-existing  turbulence  fields.  Lines  are  as 
in  Figure  1.1 
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Figure  1.6  Comparison  of  computed  and  measured  (symbols)  cross-stream  turbulence  velocities  at  the 
location  of  the  obstacle  for  turbulent  caseous  explosions  in  pre-existing  turbulence  fields.  Lines  arc  as 
in  Figure  1 . 1 
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DISCUSSION 

The  results  shown  above  illustrate  the  strong  influence  of  initial  turbulence  levels  upon  the  evolution 
of  the  explosion  kernel.  A  main  benefit  of  the  use  of  turbulence  to  enhance  the  DDT  process  is  a 
reduction  in  the  relative  sensitivity  to  the  reactivity  of  the  fuel.  The  drawback  is  that  the  turbulence 
intensities  and  scales  have  to  be  such  that  rapid  flame  propagation  is  promoted  without  undue 
quenching  due  to  turbulence  a  difficult  balance  to  obtain.  Front  a  theoretical  perspective  it  is 
important  to  point  out  that  no  attempt  to  model  such  effects  in  a  rigorous  manner  has  been  reported  :n 
the  literature  to  dale.  As  a  first  step,  the  possibility  of  performing  transient  transported  PDF 
calculations  is  assessed  in  the  section  below.  It  must,  however,  be  emphasized  that  carefully 
controlled  experimentation  is  necessary  given  the  extreme  complexities  associated  with  the  flows. 

TASK  2:  FURTHER  EVALUATION  OF  1  HE  POTENTIAL  OF  1 RANSPORTED  PDF 
APPROACHES  IN  THE  CONTEXT  OF  TIME  DEPENDENT  CALCULATIONS. 

Work  on  I  csk  2  centered  on  the  extension  of  the  high  Reynolds  number  parabolic  transported  PDF 
calculations,  reported  as  pan  of  the  earlier  contract  F6I775-00-WE054.  to  include  elliptic  time- 
dependent  cases  where  boundary  conditions  result  in  large  scale  flow  instabilities.  The  predominant 
concern  here  is  the  validation  of  the  calculation  procedure  within  the  resources  available 
Furthermore,  due  to  the  need  for  detailed  comparisons  with  measurements,  computations  have 
therefore  been  preformed  for  bluff  body  stabilized  flames  at  a  range  of  Reynolds  numbers.  The  fuel 
featured  to  date  comprises  of  a  mixture  of  naniral  gas  and  hydrogen.  The  flame  under  consideration 
here  is  the  axi-symmetric  bluff  body  stabilized  diffusion  flame  investigated  extensively  over  many 
years  by  Masri  and  co-workers  at  the  University  of  Sydney  (e.g.  Dally  et  al  1998).  The  burner 
geometry  and  layout  is  shown  below.  Time-dependent  calculation  of  the  type  presented  here  are 
exceptionally  resource  demanding.  Indeed,  the  "proof  of  concept"  calculations  shown  here  are 
probably  the  most  comprehensive  performed  to  date. 
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Figure  2.2 

Instantaneous  \eiocity  and  temperature  fields  for  the  Sydney  bluff  body  burner  flame  denoted  HM1 
obtained  through  the  current  unsteady  RANS  computation  coupled  with  a  transported  PDF  approach 
closed  at  the  joint  scalar  level  for  the  thermochemistry. 


Figure  2  J 

Averaged  velocity  and  temperature  fields  for  the  Sydney  bluff  body  burner  flame  denoted  HM1 
obtained  through  the  current  unsteady  RANS  computation  coupled  with  a  transported  PDF  approach 
closed  at  the  joint  scalar  level  for  the  thermochemistry. 

I  he  flow  ts,  unsteady  and  elliptic  and  the  equation  set  and  solution  technique  for  the  How  field  is  as 
given  in  Appendix  A  with  the  key  extension  that  a  transported  PDF  approach  with  a  scalar  space  of  16 
independent,  4  dependent  and  28  steady  state  species  is  applied  to  close  the  thermochemistry  The 
chemistry  applied  is  identical  to  that  of  the  parabolic  jet  flame  studied  earlier  by  Lindstcdt  el  cl. 
(2000),  The  same  time-dependent  elliptic  calculation  procedure  featuring  a  second  order  accurate 
TVD  scheme  on  a  staggered  grid  is  used  in  both  cases.  The  calculation  is  performed  is  a  two- 
dimensional  axisymmetric  geometry  with  a  distributed  124x109  computational  grid.  Good  temporal 
resolution  :s  maintained  with  a  Cournnt  number  not  exceeding  0.1  for  the  finite  volume  part  of  the 
computation.  The  solution  procedure  for  the  PDF  equation  features  between  200  and  400  moving 
Lagrangian  particles  per  cell.  In  addition,  a  variable  number  of  steps  are  used  to  integrate  the  PDF 
equation  during  each  finite  volume  step  and  the  mean  density  field  is  filtered  before  it  is  passed  to  the 
finite  volume  solver.  A  key  poinl  of  interest  is  the  evident  ability  of  the  approach  to  correctly  identify 
the  large  scale  instability  in  the  upper  shear  layer  observed  experimentally.  The  very  significant 
differences  between  a  realization  of  the  instantaneous  and  averaged  velocity  fields  is  readily  apparent 
through  a  comparison  of  Figures  2.2  and  2.3.  The  resulting  mean  velocity  fields,  extracted  by 
averaging  the  calculation  for  part  (~  50“  c.i  of  an  instability  cycle  time  can  be  seen  in  Figures  2  4  to 
2.7.  It  is  important  to  note  that  resource  limitations  precluded  a  full  integration  over  several  instability 
cycles  arid  the  results  shown  in  Figure  2.4  to  2. 1 1  arc  therefore  principally  used  for  illustration 
purposes  and  to  provide  comparisons  with  subsequent  calculations. 
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Figure  2.4 

Averaged  axial  velocities  for  the  Sydney  bluff  body  burner  flame  denoted  HMl  obtained  through  the 
current  unsteady  RANS  computation  coupled  with  a  transported  PDF  approach  closed  at  the  ;oim 
scalar  level  for  the  thermochemistry. 
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Figure  2.5 

Averaged  cross-stream  velocities  lur  the  Sydney  bluff  body  burner  flame  denoted  HMl  obtained 
tlirough  the  current  unsteady  RANS  computation  coupled  w  ith  a  transported  PDF  approach  closed  at 
the  joint  scalar  lev  el  for  the  thermochemistry. 
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Figure  2.6 

Axial  turbulence  velocities  for  the  Sydney  bluff  body  burner  flaine  denoted  HM1  obtained  through 
the  current  unsteady  RANS  computation  coupled  with  a  transported  PDF  approach  closed  at  the  joint 
scalar  level  for  the  thermochemistry. 
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Figure  2.7 

Cross-stream  turbulence  velocities  for  the  Sydney  bluff  body  burner  flame  denoted  HMI  obtained 
through  the  current  unsteady  RANS  computation  coupled  with  a  transported  PDF  approach  closed  at 
the  joint  scalar  level  for  the  thermochemistry'. 
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The  solution  of  the  PDF  Transport  equation  (with  a  closure  at  the  joint  scalar  level)  is  obtained,  as 
stated  above,  using  Lagrangiaii  stochastic  particles  and  the  modified  C  url's  mixing  model.  The  current 
fuel  is  C1LH;  <T.lt  and  the  jet  velocity  is  1  IS  ms  (IIMIE)  with  a  coflow  velocity  is  40  ms  I  he 
computation  is  performed  on  two  networked  Intel  Xcon  2.2  GHz  dual  processor  Dell  Precision  530 
workstations  using  message  Passing  Interface  (MPi)  for  inter-processor  communication  and  the 
Scalable  Parallel  Random  Number  Generator  (SPRNG)  for  generation  of  uncorrclatcd  random 
number  streams.  F.ach  available  processor  is  assigned  an  equal  number  of  particles  and  inter-processor 
communications  required  for  the  computation  of  the  mixing  process  and  compilation  of  global 
statistics.  The  speed-up  is  1.94  and  3.84  ov  er  two  and  four  processors  and  the  parallel  efficiency  can 
be  improved  further  by  using  a  better  network  configuration.  Examples  of  scalar  field  predictions  are 
given  in  Figures  2.8  and  2.9  below.  The  agreement  with  respect  to  the  final  products  of  combustion  is 
very  reasonable  though  averaging  over  several  instability  cycles,  as  outlined  below,  is  necessary  for 
full  quantitative  comparisons  with  measurements.  The  results  obtained  are.  however,  encouraging. 
Finally,  results  obtained  for  the  mean  and  rms  of  the  mixture  fraction  are  shown  in  Figures  2.10  and 
2.11.  The  work  performed  on  the  coupling  of  transported  PDF  methods  to  transient  calculation 
methods  establish  beyond  doubt  that  such  simulations  arc  possible.  The  application  of  this  technique 
to  the  studs  of  turbulence  enhanced  DDT  is  distinctly  possible  and  can  be  expected  to  prov  ide  a 
powerful,  albeit  computationally  expensive,  tool. 


Figure  2.8 

Mean  carbon  dioxide  concentrations  for  the  Sydney  bluff  body  burner  flame  denoted  HMl  obtained 
through  the  current  unsteady  RANS  computation  coupled  with  a  transported  PDF  approach  closed  at 
the  joint  scalar  level  for  the  thermochemistry. 


Coupled  transient  simulations  have  also  been  performed  covering  time-averaging  of  results  covering 
several  instability  cycles  to  enable  more  reliable  comparisons  with  measurements.  Computations  were 
restarted  from  the  uncoupled  calculations  by  Lmdstedt  and  Kuan  (2CmU) 
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Figure  2.9 

Mean  water  concentrations  for  the  Sydney  bluff  body  burner  flame  denoted  HMI  obtained  through 
the  current  unsteady  RANS  computation  coupled  with  a  transported  PDF  approach  closed  at  the  joint 
scalar  level  for  the  thermochemistry. 
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Figure  2.10 

Mean  mixture  fraction  profiles  for  the  Sydney  bluff  body  burner  flame  denoted  HMI  obtained  through 
the  current  unsteady  RANS  computation  coupled  with  a  transported  PDF  approach  closed  at  the  joint 
scalar  lev  el  for  the  thermochemistry. 
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Figure  2. 11 

RMS  mixture  fraction  profiles  for  the  Sydney  bluff  body  burner  flame  denoted  HM I  obtained  through 
the  current  unsteady  RAMS  computation  coupled  with  a  transported  PDf  approach  closed  at  the  joint 
scalar  level  for  the  thermochemistry. 


Time  averaging  performed  according  to  the  simple  method  by  Muradog'  i  el  ul.  (1999)  with  the 
density  evaluated  as 


'l  —  Ip"’'  +  -  —  p' 
V  V 

- »  jy  /  ‘Mi 


where  ,V7  is  the  under-relaxation  factor. 

Time  averaging  is  useful  in  reducing  the  effects  of  fluctuations  induced  by  statistical  errors  in 
the  density  field  w  ithout  either  increasing  the  number  of  stochastic  particles  or  the  use  of  statistical 
resampling.  How  ever,  the  temporal  resolution  and  the  physics  of  the  problem  is  affected.  Smoothing 
algorithms,  such  as  cross- validated  cubic  smoothing  splines  can  be  problematic  particularly  in 
evolving  regions  with  sleep  gradients  C  ross  validated  cubic  splir.e  techniques  use  statistical 
information  to  determine  the  amount  of  smoothing  required  and  thus,  arc  also  subject  to  statistical 
error.  In  other  approaches,  not  featuring  cross-validation,  the  degree  of  smoothing  needs  to  be 
specified.  Furthermore,  smoothing  does  little  to  improve  the  statistical  accuracy  of  the  data  set. 

Discrete  Monte  Carlo  sampling  converges  with  a  rate  proportional  to  the  square  root  of  the 
number  of  stochastic  panicles.  The  statistical  error  can  be  reduced  by  statistical  resampling  which  can 
be  achieved  by  increasing  the  number  of  Monte  Carlo  steps  (with  reduced  time  step  for  each  Monte 
Carlo  step)  within  each  fractional  step  of  the  coupled  calculation  procedure.  The  approach  reduces 
ihe  memory  requirements  but  docs  little  to  alleviate  the  computational  time  required. 

Computations  with  no  time  averaging  shown  above  illustrate  that  the  directly  coupled 
transient  calculation  procedure  works  well  in  principle.  A  complication  with  the  calculation  reported 
here  is  that  a  larger  length  scale  estimate  was  applied  at  the  bluff  body  boundary  This  choice  of 
boundary  condition  results  in  much  more  pronounced  vortex  shedding  which  strongly  affects  the 
upper  shear  layer  and  the  structure  of  the  recirculation  zone.  Similar  sensitivities  have  been  observed 
in  the  LES  study  of  KcmpftV  ul.  (2004). 
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('m  en  the  uncertainties  in  the  boundary-  conditions  and  statistical  resolution,  a  strongly  under- 
relaxed  t  Y  1000)  simulation  was  performed.  By  comparison,  the  results  shown  above  in  figure 
2.4  to  2  1  were  obtained  with  A,-.-  1.  The  former  computations  features  the  same  boundary 

conditions  a>  discussed  by  Kempf  tv  ai.  i2004).  Although  Muradoglu  et  a /  (1999)  used  time 
averaging  in  the  context  of  a  steady  -talc  solution  procedure,  the  approach  is  here  extended  to  a 
transient  calculation.  A  temporal  resolution  of  the  order  of  I04  s  is  arguably  required  to  resolve  the 
periodic  fluctuations  encountered  in  the  current  flow.  In  the  simulation,  the  time  steps  are  of  the  order 
of  10  s  which  results  in  a  temporal  resolution  of  the  order  of  10'  s.  Despite  the  above  reservations,  a 
transient  behaviour  is  observed  and  further  details  are  given  below .  Evidently,  it  is  desirable  to 
provide  hetter  temporal  resolution  and  to  reduce  statistical  errors.  The  two  methods  which  appear 
most  favourable  are  to  simply  increase  the  number  of  stochastic  particles  and/or  to  use  statistical 
resampling  Better  smoothing  algorithms  may  also  be  useful. 

The  evolution  of  the  various  scalars  (temperature  and  species  mass  fractions  of  Oil.  H:0,  CO, 
CO:  and  NO)  is  illustrated  through  the  time  histories  in  Figure  2  12  After  an  initial  lime  window  of 
approximately  6  ms.  regular  periodic  fluctuations  are  observed  at  locations  close  to  the  bluff  body  fx 
60  mmi  Further  downstream,  the  flame  extinguishes  due  to  large  scale  fluid  Structures  (see  Figure 
2.13)  after  ihc  restart  and  is  still  in  the  process  of  reuniting.  A  larger  integration  time  is  therefore 
required  and  the  instantaneous  contour  plot  at  10.S  ms  after  the  restart  can  be  found  in  Figure  2. 14 

The  results  presented  below  are  averaged  over  a  time  window  of  4.5  ms  beginning  6  ms  after 
the  restart  \s  can  be  seen  from  Figure  2.12.  4  periodic  cycles  with  a  period  of  approximately  1.2  ms 
(833  Hz )  can  be  observed  in  this  time  window.  The  scalar  field  predictions  arc  compared  to 
measurements  in  Figures  2.15  -  2.22  for  axial  locations  with  x  <  65  mm.  Results  for  locations  further 
downstream  arc  not  shown  as  the  flame  is  still  developing.  The  results  are  similar  to  that  of  the 
uncoupled  simulation  although  the  predictions  are  arguably  better  at  the  first  three  measuring  stations 
The  predicted  temperature  at  x  65  mm  is  lower  than  measurements  and.  as  discussed  above,  further 
temporal  integration  is  required. 
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Figure  2.12  Evolution  of  various  scalars  after  initialization  from  the  uncoupled  transported  PDF 
solution.  The  black,  red.  green  and  blue  lines  represent  time  histories  at  locations  1,  2,  3  and  4 
respectively  as  marked  in  Figure  2. 1 3  In  the  simulation  shown  .V  .  1 000 
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Figure  2.13  Contour  plot  of  the  time  averaged  temperature  [KJ  for  flame  HMI  The  locations 
correspond  to  the  time  histories  in  Figure  2  12. 
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Figure  2.14  Contour  plot  of  the  instantaneous  temperature  [K]  for  flame  HMI  at  10.8  ms  after  the 
computation  has  been  restarted  The  locations  correspond  to  the  time  histones  in  Figure  2.12. 
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Figure  2.15  Radial  profiles  of  the  lime  averaged  mean  mixture  fraction  for  flame  HMI.  The  symbols 
represent  experimental  data  and  solid  lines  predictions. 
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Figure  2.16  Radial  profiles  of  the  time  averaged  mixture  fraction  rms  for  flame  HM1.  The  symbols 
and  lines  are  as  in  Figure  2. 1 5 
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Figure  2.'-  Radial  profiles  of  the  time  averaged  temperature  for  flame  HMI  The  symbols  and  lines 
are  as  in  Figure  2.15. 


Figure  2. IS  Radial  profiles  of  the  time  averaged  OFi  mass  fraction  for  flame  HMI.  The  symbols  and 
lines  aie  as  in  Figure  2.15. 
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Figure  2.19  Radial  profiles  of  the  time  averaged  HO  mass  fraction  for  flame  HMl  The  symbols  ar.d 
lines  are  as  in  Figure  2. 1 5 


Figure  2.20  Radial  profiles  of  the  time  averaged  CO  mass  fraction  for  flame  HMl.  The  symbols  and 
lines  are  as  in  Figure  2.15. 


Figure  2.21  Radial  profiles  of  the  time  averaged  CO:  mass  fraction  for  flame  HML  The  symbols  and 
tines  arc  as  in  Figure  2.15. 
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Figure  2.22  Radial  profiles  of  the  time  averaged  NO  mass  fraction  for  flame  HM1.  1  he  symbols  and 
lines  are  as  in  Figure  2.15 

CONCLUSIONS 

The  work  reported  here,  as  part  of  the  current  contract  (F6I775  -  01-WE050)  extends  past  efforts  in 
two  directions.  The  first  is  related  to  the  sensitivity  of  the  initial  explosion  phase  l  prior  to  DDT)  to  the 
state  of  the  mixture  resulting  from  injection  of  the  reactant  mixture.  T  he  work  naturally  highlights  the 
need  to  reconcile  the  limiting  source  term  closures  controlling  the  heat  release  as  outlined  as  part  of 
previous  work.  The  second  aspect  of  the  current  work  features  computations  of  two-dimensional 
unsteady  flows  with  comprehensive  chemistry  and  a  transported  PDF  approach  closed  at  the  joinr 
scalar  level  for  the  thennochcmistrs  This  part  of  the  work  is  exceptionally  resource  intensive  and 
constitutes  an  essential  "proof  of  concept '  study  aimed  at  providing  a  computational  procedure  in 
which  the  chemistry  of  ihe  fuel  appears  naturally  in  closed  form.  The  study  does  require  significant 
computational  resources  and  the  test  case  has  therefore  been  chosen  to  be  computationally  tractable 
within  the  current  framework,  The  results  obtained  illustrate  that  the  approach  is  technically  possible 
and  constitutes  a  potential  route  for  dealing  with  the  chemical  kinetic  effects  associated  with 
detonation  initiation  in  a  Strongly  turbulent  flow  field.  However,  it  must  be  recognised  that  more 
extensive  studies  are  required.  Permission  to  compute  the  actual  geometry  via  IT  A  R, 'Export 
Controlled  information  permission  had  to  be  abandoned  in  favour  of  generic  calculations  exploring 
the  influence  of  turbulence  levels  on  the  strength  and  timing  of  the  onset  of  the  DDT  process  The 
results  show  a  strong  influence  and  indicate  a  possible  solution  to  the  problem  of  detonation  initiation. 
It  is.  however,  also  recognised  that  computations  of  the  sort  attempted  stretch  our  current 
understanding  of  turbulent  combustion  to  the  limit.  Accordingly,  experimental  studies  that  can  further 
highlight  the  physics  are  strongly  recommended. 

DEVIATIONS  FROM  WORKPLAN 

[  he  initial  delay  caused  by  staffing  difficulties  with  respect  to  students  appeared  to  be  resolved 
satisfactorily  until  April  2002  when  the  student.  Mr.  <j.  Colangelo.  working  on  the  project  had  to 
return  to  Italy  for  health  reasons.  Further  delays  were  caused  by  the  difficulties  in  obtaining  an  export 
license  application  is  being  pursued  by  the  Air  Force  via  the  EOARD  in  order  to  advance  the  final 
element  of  Task  1  -  which  entailed  the  computation  of  geometries  of  direct  relevance  to  the  Air  Force. 
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ABSTRACT  1  he  transition  of  an  initially  laminar  or  turbulent  deflagration  to  a  confuted  gaseous  explosion  or 
detonation  i  DDT i  is  of  fundamental  relevance  to  practical  devices  such  as  pulsed  detonation  engines  I  PDF's  t  and 
remains  of  key  importance  to  on-  and  off-shore  hazard  assessment  procedures  "he  topic  is  also  one  of  the  most 
challenging  ir.  contemporary  physics  and  requires  a  detailed  consideration  of  the  interactions  between  unsteady 
high  speed  turbulent  flows  and  chemical  kmet.es  Furthermore,  in  many  cases  fluid-wall  interactions  come  to  the 
fore.  It  is  therefore  perhaps  not  surprising  that  the  subject  has  traditionally  been  treated  cither  on  the  basis  of 
shock  dynamics  n  laminar  flows,  typsea.  y  combined  with  simplified  chemica.  kinetic  expressions,  or  as  a 
turbulent  reacting  flow  closed  at  the  eddy  viscosity  level  in  which  Ihe  chemistry  :s  considered  to  be  in  the  high 
Damkohler  number  limit.  The  limitations  associated  with  each  approach  are  obvious  and  experimental  data  sets 
assembled  over  many  decades  show  -  without  doubt  -  that  the  combined  influence  of  chemistry  and  flow  mutter-, 
to  the  DDT  process  and  by  implication  neither  treatment  is  satisfactory  The  -.  urrent  chasm  in  the  theoretical 
treatment  of  turbulent  gaseous  explosions  and  detonation  propagation  studies  must  be  bridged  in  Order  to  provide 
a  technically  meaningful  modelling  capability  for  turbulence  enhanced  DDT  processes.  The  focus  of  the  present 
paprr  is  to  use  detailed  time- resolved  experimental  data  sets  to  explore  the  ability  of  comprehensive  moment 
based  closures  to  reproduce  time  resolved  flow  field  features  up  to  the  point  of  onset  of  detonation.  A  new 
reaction  rate  closure  ts  also  derived  and  evaluated  The  experimental  dai.-.  set  considered  features  a  confuted 
obstacle  accelerated  premixed  turbulent  flame  of  physical  relevance  to  PDFs  Tune-resolved  mean  and  rms 
profiles,  obtained  using  laser  Doppler  anemometry,  arc  used  for  comparisons  with  computational  results  along 
w  ith  pressure  traces  and  instantaneous  spark  -chlieren  photographs.  The  main  conclusions  of  the  study  are 
positive  and  the  work  does  show  that  the  application  of  higher  moment  based  cioi-ures  to  model  the  flow  field 
leading  up  to  the  onset  of  DDT  yields  surprisingly  satisfactory1  results 
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IM  KODl  C  TION 

Calculation  methods  tor  transient  turbulent  reacting  flows  must  balance  the  fidelity  of  thermoebeinieal  and  fluid 
mechanical  closure  aspects  in  order  to  provide  technically  useful  results  The  confined  htgh  speed  combusting 
Hows  considered  here  are  dominated  by  flow  sr.icle  interactions  and  present  a  canicular  challenge,  Chemical 
reaction  occurs  in  the  fine  (unresolved  I  scales  and  must  be  modelled  irrcspccto  e  of  the  technique  adopted  for  the 
flow  field  Furthermore,  flow  interactions  with  solid  boundaries,  c.g.  in  the  context  of  shock  (turbulent 
boundary  layer  interactions,  is  a  severe  challenge  for  both  unsteady  Reynolds  Averaged  Xavier- Stokes 
approaches  IRAN'S  and  Large  F.ddy  Simulations  iLK.S)  even  in  the  absence  ?f  chemical  reaction  I emporal 
response  issues  and  flovv-field-Turbulence  interactions  with  chemistry  rose  occupy  central  rotes  during  the  DDT 
phase  in  pulsed  delegation  engines  [1.2]  and  in  emerging  propulsion  devices  m  general  [3],  The  former  arc.  In 
contrast  lo  the  majority  of  contemporary  propulsion  dev  ices,  intrinsically  unsteady  and  cover  an  unusually  wide 
spectrum  of  thermochemistry'  and  physics  hxamples  of'  topics  of  practical  relevance  include  fuel  mixture 
injection  and  obstacle  (or  shock)  enhanced  localized  turbulent  explosions  resulting  in  a  transition  to  detonation 
followed  by  the  propagation  of  the  resulting  detonation  and  its  interaction  with:  confinement  boundaries 

Key  challenges  that  remain  to  be  resolved  during  the  latter  stages  of  the  DDT  process  include  shock  interactions 
with  flames  and  boundary  layers.  However,  two  of  the  principal  practical  issues  at  present  arguably  concern  the 
relative  influences  of  fuel  structures  and  turbulence  on  the  transition  process  leading  up  to  the  actual  onset  of 
detonation  The  chemical  structure  of  the  fuci  h.»  been  shown  experimentally  to  have  a  leading  order  influence 
on  the  transition  process  in  smooth  tubes  [4  It  ts  also  well  established  that  rite  influence  of  turbulence  reduces 
the  scnsitiv  uy  of  the  transition  process  to  the  point  where  alkane  fuels,  wnh  trie  exception  of  methane,  show  a 
similar  sensitivity.  Lower  alkenes  and  alkynes  do,  however,  retain  a  greater  propensity  to  DDT  even  in  strongly 
turbulent  environments  [5],  An  exceptionally  wide  range  of  turbulence  Reynolds  numbers  is  normally 
encountered  during  the  transition  process  and  the  resulting  complex  interactions  between  flow  and  chemistry  j-c 
likely  to  play  3  key  role  in  device  design.  Ihc  resulting  demands  placed  on  corrptr.-.uonal  design  techniques  arc 
severe.  Similar  issues  have  for  a  long  lime  prevailed  in  the  context  of  hazards  related  studies  of  relevance  to  the 
modelling  of  turbulent  gaseous  explosions.  Much  of  the  work  in  the  latter  area  h3s  been  performed  in  the  context 
of  successive  dies  sponsored  by  the  CEC  iDGXIli  as  pan  of  the  Major  1  echnological  Hazards  Programme 
and  a  summary  of  some  of  the  results  has  been  presented  by  Amtzen  t-t  at  |6J. 

The  modelling  of  turbulent  gaseous  explosions  has  to  date  exclusively  featured  different  variants  of  closures  at 
the  eddy  viscosity  level  Generally,  such  closures  have  only  proved  adeqjaic  when  turbulence  ts  driven  by  K  ca 
cross-stream  gradients.  In  turbulent  reacting  flows,  by  contrast,  variable  density  effects  are  at  least  comparable 
to.  and  usually  outweigh,  such  generation  mechanisms.  A  consequence  of  chemical  reaction  ts  the  occurrence  of 
velocity  divergence  leading  to  a  reduction  of  turbulence  levels  due  to  the  enhanced  contribution  of  the  strain 
components  m  the  normal  Reynolds  stress  “preducuon"  terms.  However.  \i  ame  expansion  will  also  result  in 
production  through  preferential  acceleration  effects  through  mean  pressure  gradient  terms  in  the  respective 
transport  equations.  Similar  considerations  apply  to  the  flux  vector  which  may  not  necessarily  be  aligned  with 
that  of  the  scalar  gradient  and  therefore  result  in  ''non-gradient"  transport.  Effects  of  the  type  outlined  above  can 
not  he  me.  rp< rated  into  eddy  viscosity  closures  Advantages  of  higher  moment  closures  are  that  “history"  and 
convection  related  effects  on  the  evolution  of  pjrbuler.ee  are  incorporated  in  a  natural  manner,  while  production 
terms  by  mean  strain  components  appear  in  a  closed  form.  Furthermore.  effects  induced  by  anisotropy 
generating  mechanisms,  such  us  mean  pressure  gradients,  are  also  included.  Many  of  die  latter  features  can  be 
expected  to  be  critical  for  the  flows  considered  here  as  it  may  be  expected  that  pressure  gradient  effects  will  have 
a  significant  impact  upon  the  evolution  of  turbulence.  I  he  application  of  higher  moment  based  closures  to  flow  s 
with  strong  pressure  gradients  is,  however,  not  trivial  and  presents  significant  numerical  challenges 

In  the  current  context  it  can  also  be  expected  that  techniques  capable  of  dealing  with  direct  kinetic  effects  -  with 
large  deviations  from  equilibrium  -  will  be  required  in  order  to  model  the  full  DDT  process.  At  present,  the 
transporter  PDF  approach  is  the  only  method  by  which  this  can  be  achieved.  Hulek  &  Lindstedt  f 7]  have  shown 
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that  one-dimensional  transient  premixed  turbulent  flames  can  be  computed  with  a  transported  PDF  closure  at  the 
joint  velocity-scalar  level.  The  molecular  mixing  terms  was  dosed  using  a  binomial  Langevm  model  modified 
for  joint  velocity-scalar  statistics  [8  ].  Lmdstedt  Ac  Vaos  |0  subsequently  exp  ore  J  die  ctTecLs  of  different  nrntic 
models  on  the  computed  turbulent  burning  ve'.> city  and  showed  that  both  the  binomial  I  angevin  and  Cur  \ 
mixing  mode  s  perform  well  when  combined  with  a  closure  at  the  second  moment  level  In  both  cases,  the 
chemical  reaction  rate  source  term  was  extracted  from  laminar  flame  calculations  and  reduced  to  a  one-step 
formulation  Time  dependent  calculations  with  realistic  (large)  scalar  spaces  (~  2(>  independent  scalarsi  are 
becoming  increasingly  feasible  and  are  starting  to  appear  [10]  though  exceptional  computational  resources  are 
required 

In  the  present  work  high  Damk&hlet  number  techniques  are  applied  to  explore  the  ability  of  higher  mon  ent 
based  closures  for  the  flow  field  to  reproduce  major  flow  field  fea  tres  i  g  lip  to  the  point  of 
detonation.  Cathn  &  I  mdsteilt  1 1 1 1  have  shown  that  reaction  rate  expressions  based  on  such  assumptions  arc 
prone  to  numerical  inabilities  at  the  leading  edge  of  turbulent  premixed  flames  I  he  result  is  that  the  genera  iy 
accepted  leading  order  scaling  between  the  turbulent  burning  velocity  and  turbulent  velocity  fluctuations  is  no; 
reproduced  m  numerical  calculations  unless  a  Heaviside  function  is  introduced  2t  the  leading  edge.  A  further  key 
observation  s  that  the  burning  velocity  eigenvalue  is  strongly  dependent  upon  such  modifications  and  that 
common  approaches  can  he  readily  expected  to  lead  to  errors  in  excess  >  t  100  \  further  observation  made  i> 

that  the  turbulent  flame  brush  thickness  has  :o  be  well  resolved  computationally  in  order  to  ensure  physically 
meaningful  results.  The  adherence  to  such,  comparatively  simple,  requirements  ensures  that  the  intrinsic  scaling 
behaviour  of  the  solved  equations  is  retained.  A  recent  study  by  Cotangeio  A  i.indstedt  ( 1 2]  has  shown  that  the 
temporal  evolution  of  gaseous  explosions  can  be  captured  with  usefu.  accuracy  ar.d  that  alternative  reaction  rate 
closures  based-  on  3  fractal  assumption  with  the  limiting  cases  of  inner  cut-off  scales  corresponding  to  the 
Kolmogorov  ar.d  Gibson  scales  bracket  key  aspects  of  the  experimental  data 

Reliable  experimental  data  is  exceptionally  difficult  to  procure  for  strongly  transient  turbulent  flows  but  remain  a 
prerequisite  for  a  scientific  evaluation  of  different  closure  elements  The  standard  manner  of  evaluating 
predictive  techniques  with  reference  to  single  parameters,  such  a  pressure  rise  cr  propagation  velocity,  is  3iso 
insufficient  Furthermore,  while  visual  information  is  helpful  in  advancing  qualitative  understanding,  any 
detailed  assessment  of  controlling  mechanisms  requires  quantitative  flow  field  information.  It  is  therefore  a 
striking  deficiency  that  at  the  present  lime  only  two  sueli  data  sets  exist  at  meaningful  Reynolds  number'. 

1 1 3.14,1 5.16]  The  current  work  takes  advantage  of  one  of  these  data  sets  to  provide  an  assessment  of  the  ability 
of  high  Pamkf'hicr  number  based  approaches  to  model  qualitatively  and  quantitatively  the  evolution  of  a 
turbulent  gaseous  explosion  tn  a  confined  channel  The  case  considered  features  a  pre-existing  turbulence  field 
generated  through  premixed  fuel-air  mixture  injection  at  the  closed  end  or'  a  detonation  tube  [15,16  ;.  In  the 
present  work  closures  at  the  second  moment  level  are  util i7ed  for  both,  velocity  and  scalar  fields  The 
computational  techniques  applied  are  capable  of  describing  phenomena  such  as  flow  field  anisotropy  and  o: 
providing  qualitatively  correct  flame  dynamics  leg  Lindstcdt  &  Vaos  [17]}  The  recent  developments  of 
comprehensive  closures  at  this  level  [l'.l*]  is  here  followed  by  their  application  to  model  u  transient  multi¬ 
dimensional  premixed  turbulent  flame 
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(.0\  LRMM.  EOLATIONS 

The  aerothcnm-chemieal  fields  are  governed  by  the  conservation  equations  tor  muss.  momentum  anti  scalars  cast 
in  a  density  weighted  form  I  aminar  diffusion  terms  .aid  extenial  force  fields  are  neglected. 
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In  the  context  of  time  dependent  flows  it  is  important  to  note  that  changes  in  the  velocity  and  scalar  fieids  drives 
die  generation  of  turbulence  and  that  die  feedback  is  provided  though  the  turbulent  correlation  terms  on  the  RHS 
of  equations  (2.3 1  The  latter  are  here  treated  through  the  corresponding  transport  equations  for  the  Reynolds 
Stresses  (4-6)  and  scalar  (luxes  (7-8). 
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I  he  terms  on  The  R1IS  represent,  in  order,  turbulent  transport  of  the  Rcyr.o.ds  stresses,  effects  of  mean  strain  tor 
"production"  P. 1  effects  of  mean  pressure  gradients  the  turbuleni  pressure  strain  term  (fc.)  and  viscous 
dissipation  ipt.  > 
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The  modelling  of  the  pressure  correlation  terms  n  the  Reynolds  stress  equations  is  the  focal  point  ot  second 
moment  closures  The  terms  am  be  decomposed  into  redistributive  ami  isotropic  parts  and  details  can  be  found 
elsewhere  ic.g.  Lindstcdt  &  Vaos  fl7;l  The  i-.tter  may  be  further  split  into  pressure  transport  and  pressure 
dilution  terms.  Redistribution  terms  are  invariably  modelled  by  recasting  closure's  (derived  on  a  constant  density 
basis)  in  a  density  weighted  form  and  a  sim.iar  procedure  is  also  applied  to  the  pressure  scrambling  terms.  An 
algebraic  expression  for  rhe  density  weighted  fluctuations  forming  part  of  the  mean  pressure  gradient  term  can 
be  dcriv cd  T he  corresponding  scalar  (luxes  for  the  reaction  progress  variable  are  given  below. 
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I  he  present  study  considers  fully  compressible  (lows  with  compressibility  effects  introduced  via  the  solution  ot 
a  conservation  equation  for  the  internal  energy  and  the  equation  of  state.  The  thermodynamic  properties  are 
evaluated  from  JaNaF  polynomials  and  the  temperature  is  calculated  via  a  Newton  iteration. 


(9)  ape  0 pu,  e  _r»w,  dpu,e  -  7.  u  a? 

— — + - -  - -  pu.e  - and  PM  =  pRT  i  -  i 

or  ax.  Ox,  Os  <7,  Ox, 


The  treatment  ot  flux  terms  in  the  interna  energy  equation  presents  a  currently  unresolved  challenge  and  the 
principal  roic  of  Eq.  (9)  in  the  present  context  is  to  provide  a  treatment  for  the  calculation  ot  the  temperature  and 
density  field  I  he  turbulent  kinetic  energy,  required  n  the  evaluation  of  the  “eddy  viscosity”,  is  determined  from 
the  solution  .  ■*'  the  Reynolds  stress  equations  and  the  turbulent  I’randt)  number  is  assigned  a  value  of  0.  ;5. 
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The  treatment  ..f  the  pressure  strain  term'  features  a  rearrangement  into  redistributive  and  isotropic  parts  with  the 
latter  decomposed  into  pressure  transport  and  dictation  terms. 
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The  corresponding  pressure  strain  term  in  the  scalar  (lux  equation  may  he  re.j't  in  terms  of  pressure  transp  rt 
t T>  and  scrambling  IS)  terms. 
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The  specific  mode:  expressions  applied  are  described  below. 


specific:  modei  forms 

Tlie  “return  to  isotropy"  anti  “strain"  redistribution  par.s  are  obtained  from  model  formulations  derived  in  non- 
reacting  flows.  Several  suggestions  exist  a:  various  levels  of  complexity.  The  model  used  here  for  illustrative 
purposes  is  that  of  Haworth  &  Pope  [19]  which  is  based  on  a  stochastic  Lagrangian  Generalized  I  ungevir,  Model 
(GLM)  [20]  The  latter  features  a  non-linear  return  to  isotropy  and  has  been  calibrated  and  tested  for  a  range  •  f 
homogeneous  turbulence  and  free  shear  flow*  An  important  feature  of  the  model  is  That  :t  also  yields  a  closure 
for  the  pressure  scrambling  term  in  the  scalar  flux  equation  and  thus  allows  for  realizable  modelling  of  both 
pressure  strain  correlations.  The  "return"  and  ''strain"  redistribution  parts  take  the  following  equivalent  forms 
where  b  denotes  the  anisotropy  tensor. 
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The  corresponding  model  form  for  the  scalar  flux  equation  is 
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In  Eq  (13 1  the  constant  Cu  is  assigned  the  value  2.1  fhe  functional  form  tor  (/  <  linear  in  the  mean  velocity 

gradients  and  quasi-lmear  in  the  anisotropy  tensor  and  the  expanded  form  :s  given  in  Appendix  A.  A  number  of 
closures  have  been  suggested  for  the  triple  moments.  However,  the  contributions  of  the  diffusion  terms  to  the 
overall  budget  '  the  Reynolds  Stresses  and  scalar  fluxes  is  moderate  and.  consequently,  the  turbulent  transport 
of  the  second  moments  is  here  modelled  using  the  generalized  gradient  diffusion  traxlel  of  Daly  &  1  lut  low  [21. 
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The  closure  for  the  dissipation  rate  of  the  turbulent  kinetic  energy  is  obtained  by  the  standard  variable  density 
transport  equation  model  (e.g.  [17,22]). 


CLOSURE  FOR  REACTION  RATE  REL  \TED  TERMS 

The  closure  for  the  chemical  source  term  presents  an  interesting  problem  that  entails  the  determination  of  an 
appropriate  time  scale  for  scalar  turbulence  Traditionally,  the  closure  for  the  scalar  dissipation  rate  in  the  high 
Damkohler  number  regime  is  purely  based  on  a  constant  timescale  ratio  leading  to  a  direct  relation  to  die 
mechanical  time  scale.  The  analysis  is  here  pursued  from  the  basis  of  differentia :  modes  and  the  form  proposed 
by  Jones  ic  Musenge  (23)  is  used  here  for  illustrative  purposes. 
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In  the  above  equation.  <.  ndicates  a  model  constant. ./  a  turbulent  transjvirl  term  and  the  production  terms  <P  P.  < 
follow  from  the  equations  tor  the  evolution  or  the  Reynolds  stresses  and  sca...r  fluxes  By  virtue  of  Eq  1 1“>  an 
algebraic  expression  for  the  scalar  timescale  can  be  derived  if  it  is  assumed  that  the  evolution  of  the  scalar 
dissipation  rate  s  dominated  by  local  source  terms  The  resulting  form  reduces  :o 
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The  above  procedure  may  be  regarded  as  typical  for  the  derivation  of  the  commonly  adopted  scalar  time-scale 
ratio  The  obvious  requirement  is  that  ihc  sum  of  the  terms  in  the  square  bracket  does  not  vary'  significantly 
Perhaps  surprisingly,  the  value  for  the  time  scale  ratio  ;s  typically  of  order  ur. its  for  passive  scalars  in  many 
flows  In  cases  featuring  reactive  scalars  a  similarly  derived  equatioi  for  the  scalar  dissipation  rate  includes 
additional  specif  c  terms  in  the  form  of  correlations  between  scalar  and  reaction  rate  gradients. 
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Borghi  Ac  Mantel  [24 1  have  performed  a  dimensional  analysis  and  show  that  the  above  term  is  proportional  to  the 
Dtimkohlcr  number  It  has  also  been  shown  by  Lindstcdt  &  Vaos  T 17. 1  S'],  based  on  a  fractal  analysis,  that  in  the 
latter  regime  the  tune  scale  ratio  needs  to  be  modified  to  account  for  flame  propagation.  Bray  ,-t  at.  |25  have 
recently  evaluated  a  number  of  reaction  rate  clo-ures  in  die  context  of  the  opposed  je:  geometry.  The  conclusions 
are  interesting  arid  suggest  that  the  majority  of  reaction  rare  closures  fail  to  predict  qualitative  trends  correctly 
furthermore,  u  lias  been  shown  that  irrespective  of  the  level  of  closure  adopted  -  algebraic  or  transport  equation 
based  -  current  mc-deis  could  generally  no:  be  calibrated  to  perform  satisfactorily  irrespective  of  the  choice  of 
modelling  constant! -s>  An  exception  was  found  to  be  the  fractal  based  closure  derived  by  Lindstcdt  Ac 
Sakthitaran  [2b]  and  discussed  further  by  Lindstcdt  At  Vaos  [17].  It  must,  however,  be  recognized  that  Bray  et  al. 
[25]  also  in  his  vase  show  that  there  is  a  significant  uncertainty  relating  to  the  appropriate  modelling  constant 
and  that  the  value  proposed  [17]  appeal  >  too  low.  Based  on  the  above  considerations  the  resulting  modified 
expression  for  the  time  scale  ratio  is  proposed. 
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The  reaction  rate  i.VcAO  is  the  integral  across  tl e  corresponding  | laminar ■  Dame  sheet  and  Ihc  viscosity  i.n  and 
I’randtl  number  in  i  refer  to  the  corresp:  :  .  It  properties  evaluated  U  some  inner  reference  plane  It  is 

important  to  note  that  the  above  expression  .nip  ies  that  burning  velocity  reductions  at  high  Reynolds  numbers 
arc  due  to  partial  flame  quenching  and  not  to  flame  geometry  related  considerations.  The  limitations  of  the  above 
derivation  procedure  and  the  resulting  expression  arc  obvious  but  to  a  large  degree  identical  to  those  present  in 
standard  modelling  approaches  From  a  practical  perspective,  however,  the  derived  form  ensures  a  consistent 
sealing  behaviour  for  turbulent  burning  velocities  i:  the  high  Dainkohler  number  regime  of  combustion.  For  he 
present  case,  the  constant  (?\  has  been  assigned  the  standard  value  of  I  while  the  parameter  C'\  has  been 
calibrated  to  the  value  ot  1.2  on  the  basts  of  flames  stabilized  in  an  opposed  jet  geometry  The  corresponding 
reaction  rate  expression  follows  from  the  assumpfon  of  a  bimodal  PDF. 
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The  reaction  rate  constant  C#  is  assigned  the  value  of  3.5  [ 27]  and  the  above  expression  arguably  constitutes  an 
upper  iir.it  The  velocity  reaction  rate  corrciaticr.  term  is  here  closed  through  the  assumption  of  a  bimodal  PDF 
[17], 
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Different  models  for  the  pressure  dilation  icon  have  been  proposed  by  Lindstedt  &  Hulek  [.IS]  and  /.hang  & 
Rutland  1-9;  The  former  suggestion  has  the  virtue  of  being  consistent  with  the  assumed  therm  ochcmical 
conditions  and-  while  the  latter  is  heuristic  in  nature,  ir  does  not  present  problems  n  terms  of  realizability  .  In  the 
present  work  the  term  has  been  omitted  as  the  influence  is  not  expected  to  be  large  [US 

EXPERIMENT  \L  AM)  COMPUTATION  VL  CONFIGURATION 

The  transition  of  a  turbulent  flame  to  a  gaseous  explosion  lit  a  confined  pre-existing  turbulent  flow  field  b 
co:iv:«!ered  tr.  -Jte  context  of  the  data  set  produced  by  Lindstedt  &  McCann  [15, 16].  The  configuration  yields 
over-pressures  •  irouod  200  kPa  and  peat  ilk  ekdtitt  around  40  1  ::is  both  significantly  higher  than  those 
reported  o>  1  mdstedt  &  Sakthitharan  [1  ?.!•*]  for  initially  quiescent  mixture*  -  and  the  turbulence  velocities 
reach  values  of  the  order  40  nv's.  Deflagration  to  detonation  transition  occurs  ir.  the  same  configuration  for  mote 
reactive  mixtures  than  the  stoichiometric  methane-air  case  reported  here  The  detonation  tube  features  six 
interchangeable  rectangular  sections  (72  mm  x  34  mm)  of  length  1.825  m  and  two  window  sections  of  length 
0  m.  Initial  turbulence  is  generated  bv  an  axi-symmctric  jet  of  5  mm  internal  diameter  centred  in  a  plate  a: 
the  ignition  end  of  the  flame  tube.  The  flow  field  created  by  the  injection  of  the  premixed  lucl-air  mixture  has 
been  meastireu  [  151  and  the  data  is  here  used  in  order  to  form  the  initial  se:  of  condition*;.  Ignition  is  obtained  by 
two  opposing  electrodes  placed  symmetrically  either  side  of  the  jet  centre  line.  A  5  mm  thick  obstacle  with  a 
height  of  36  mm  is  placed  downstream  of  the  jet  extt  a;  a  distance  of  415  mm  The  closure  outlined  above  is  here 
applied  to  the  above  geometry  using  a  locally  refined  two-dimensional  computational  grid  (-  1()0  000  nodes) 
stiflictent  to  provide  a  grid  independent  solution.  The  computations  feature  a  second  order  TVD  scheme  and  an 
implicit  predictor-corrector  method  with  time  splitting  error  control  [11,17]. 

Ihe  applied  diagnostics  include  sixteen  coaxial  ionisation  probes,  used  to  detect  flame  arrival  times,  and  six 
piezoelectric  pressure  transducers  pos.t.oned  along  the  channel.  The  flame  propagation  process  was  aiso 
visualised  by  means  of  spark  schlieren  photography  featuring  a  parallel  beam  of  white  light.  Velocities  were 
obtained  using  the  standard  dual-beam  forward  scatter  laser  Doppler  anemometer  arrangement  and  the  dau 
reduction  technique  reported  by  I  indstedt  &  Sakthitharan  [13]  applied.  Velocity  measurements  along  the  axial 
(u  -velocity)  ar.d  vertical  directions  (v  velocity )  have  been  obtained  at  33  points  above  and  downstream  of  the 
obstacle.  All  the  points  lie  on  the  vertical  plar.e  passing  through  the  axis  of  the  flame  tube.  Sakthitharan  1 14  has 
shown,  through  flow  field  measurements  across  the  detonation  tube,  that  the  flow  is  predominantly  nvo- 
dimensiona.  in  the  region  close  to  the  obstacle  A'lhough  there  arc  temporal  variations  associated  with  the  initial 
kernel  growth  (-  10  %).  the  reproducibility  of  the  flame  arrival  at  subsequent  pons  is  excellent  with 
comparatively  small  variations  (-  3  %)  The  iatter  indicates  that  ihe  development  of  the  vortex  structure 
downstream  or  the  baffle  is  highly  reproducible  and  relatively  insensitive  to  initial  variations  Problems  do 
naturally  arise  with  measurements  of  mean  and,  in  particular,  turbulence  velocities  m  the  current  strongly 
accelerated  flow  It  :s  therefore  expected  that  the-  measurement  accuracy  is  significantly  reduced  following  the 
onset  of  the  explosion 

RESU  LTS  AND  DISCUSSION 

Computed  flame  contours  are  compared  with  schlieren  images  in  Figs  I  and  2.  The  relative  times  in  the  flame 
propagation  cycle  are  close  but  not  matched  exactly  and.  furthermore,  the  SCtllicren  images  arc  •Ime-nf-siehf 
across  the  flame  tube.  Caution  ts  therefore  required  in  interpreting  the  images.  However,  it  is  seen  that  the 
creation  of  large  scale  flame  folds  and  the  flame  wrapping  into  the  vortex  3re  qualitatively  reproduced  by  the 
computation  Clearly  the  flow-  features  large  scale  instabilities  along  with  strong  turbulence  generation  in  the 
shear  layer  and  re -circulation  zone  behind  the  obstacle.  The  large  scale  instabilities  ure  ot  order  10  mm  with  the 
wrinkling  of  the  flame  appearing  to  be  of  order  1  mm.  The  predicted  pressure  traces  are  compared  with 
measurements  in  Fig.  3a-b.  Close  to  the  ignition  [mint  the  agreement  is  very  satisfactory  as  shown  in  Fig  3a 
However,  i  is  also  evident  that  peak  pressures  are  somewhat  under- predicted  closer  to  the  obstacle  (Fig  3b). 
though  die  shape  of  the  pressure  trace  profile  remains  comparativ  elv  well  reproduced. 
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Further  evidtact  that  ihc  predictions  arc  qualitatively  good  and  quantitatively  reasonable  can  he  found  in  the 
axial  an.:  or.  ■"-stream  mean  velocity  profiles  shown  in  Figs.  4  and  5  for  the  p.ane  where  the  obstacle  s  located 
I  he  interact '''c  of  strong  flow  gradients  in  the  \ena  .  oHtracra  over  the  battle  with  the  explosion  generated 
pressure  pulse;  are  exceptionally  till  to  measure  and  it  is  likely  that  the  experimental  data  underestimates 
the  peaks  in  the  bulk  velocities.  However,  with  the  exception  of  the  latter,  the  profile  shapes  and  magnitudes  arc 
well  reproduced,  tor  the  axial  (u)  velocity  component.  The  computations  show  strong  oscillations  in  the  bum;  gas 
due  to  travelling  pressure  waves  (evident  in  the  pressure  traces  shown  m  Figs  3)  which  clearly  could  not  F-e 
captured  by  ...  ....asuremcnls.  The  ...  ....  the  ctos&etriean  (\  ■  component  is  much  less  iatisfu 
Howevet.  it  can  be  noted  that  the  How  velocities  range  from  well  in  execs*  nf  i  n  m  s  close  to  ihe  obstacle  it? 
1 42  mmi  to  around  10  m  s  closer  1 10  mmi  to  the  top  wall  It  may  further  be  nonced  that  mean  velocities  are  well 
reproduced  until  the  onset  of  the  explosion  Further  away  from  the  obstacle,  at  the  515  mm  plane,  the  agreement 
between  measurements  and  computations  i;  generally  satisfactory  for  the  axial  velocity  component  as  show  n  in 
Fig.  6.  However  it  may  be  noted  that  the  strongest  How  accelerations  again  present  what  probably  .> 
measurement  difficulties.  Hie  qualitative  and  quantitative  agreement  for  the  cross-stream  velocity  component  is 
very  satisfactory  as  shown  in  Fig  7  -  despite  tFic  fact  that  the  time  evolution  of  the  flow  is  extremely  complex  It 
may  also  be  noted  that  the  flow  acceleration  is  less  severe  in  the  cross-sireur.,  direction  and  lienee  greater 
confidence  may  be  placed  in  experimental  data. 

I  lie  difficuit.es  outlined  above  in  the  context  of  the  determination  of  the  bu  k  velocity  field  are  naturally  also 
present  in  the  measurements  of  the  lurbulencc  intensities.  It  is  important  to  note  that  the  vortex  shedding,  and 
hulk  flow  instabilities  do  nol  satisfy  the  common  criteria  (e.g.  irregularity i  as;octatcd  with  turbulence  (e.g 
Tcnnckcs  &  Lumlcy  (30 ) >.  A  scale  separation  associated  with  bulk  and  turbulent  motion  is  evidently  present  in 
the  current  flow  but  becomes  awkward  to  analyse  \.th  computationally  and  experimentally  The  combination  of 
a  time  window  of  500  pc.  used  experimental  y  in  order  to  gather  sufficient  statistics,  along  with  the  difficulties 
experienced  by  seeding  particles  in  following  flow',  renders  the  measurements  difficult  to  interpret.  Nevertheless, 
up  to  the  point  of  the  onset  of  the  explosion  I  16-17  ms)  there  is  fair  agreement  between  computed  and 
measures  values  of  the  axial  component  as  shown  it  Fig.  S.  Following  tie  uiset  of  the  explosion  the  axial 
turbulence  velocities  -  in  particular  -  are  very  strongly  affected  and  it  is  no-  possible  to  advance  conclusions 
regarding  the  agreement  between  computed  anv  measured  values  The  cross-stream  velocity  component  is  again 
much  less  affected  and  remains  surprisingly  well  reproduced  both  qualitatively  and  quantitatively  a;  shown  tn 
Fig.  9  A  comparison  of  the  results  shown  in  Figs  7  and  9  is  interesting  and  suggests  that  despite  the  difficulties 
associated  with  both  measurements  and  computations  the  overall  agreement  is  encouraging  It  is  notable  that  the 
rapid  changes  in  the  velocity  field  and  the  resulting  turbulence  generation  are  both  well  captured  Perhaps 
surprisingly,  the  results  indicate  that  the  computational  results  arc  iu  many  cases  close  to  experimental 
uncertainties  for  the  current  exceptionally  complex  flow 

CONCLUSIONS 

The  ability  :o  model  the  formation  and  evolution  of  DO  I  kernels  is  of  fundamental  importance  to  a  range  i 
practical  applications  >uch  as  pulsed  detonation  engines.  In  order  to  advance  the  subject  it  is  necessary  to  cxr  ore 
the  sensitivity  of  the  overall  phenomena  to  the  maior  closure  elements  of  physical  relevance  to  the  problem 
Experimental  work  has  shown  that  in  many  practical  areas  the  driving  force  for  the  transition  can  be  found  ir.  the 
interaction  of  turbulence  and  chemistry,  while  in  other  cases  it  will  reside  with  shock  (turbulent)  boundary  layer 
interactions  Both  topics  present  similar  -  if  not  close  to  identical  -  modelling  difficulties  irrespective  of  the 
technique  used  for  the  flow  field  (e.g  LES  or  unsteady  RANS).  Nevertheless,  it  is  naturally  essential  that  a 
sensible  reproduction  of  the  flow  field  is  achieved  by  the  computational  method  of  choice.  In  the  absence  of 
experimental  ilata  such  assessments  arc  conjectural  and  the  present  work  hs-  used  comparatively  detailed  data 
for  turbulent  gaseous  explosion  as  a  bast'  for  comparisons.  The  experimental  data  set  ts  naturally  subject  to 
significant  unccnairncs  resulting  from  exceptionally  difficult  measurement  conditions  resulting  from  repeated 
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flow  acceleration  deceleration  cycles  Despite  such  difficulties  the  work  has  shown  that  the  application  of 
comprehensive  ici  er  moment  closures  to  the  modelling  of  the  initial  onset  <•'  i.'DT  is  technically  possible  and. 
furthermore,  rcs.dts  in  encouraging  quantitative  and  qualitative  agreement  with  measurements,  It  may  be  of 
interest  to  note  that  good  agreement  for  other  unsteady  flows  have  also  been  obtained  le  g.  Durbin  [ 33] I  using 
moment  based  methods.  Ihe  turbulence  chemistry  interactions  can  in  principle  be  closed  using  a  transported 
PDF  approach  combined  with  comprehensive  autc-igntnon  chemistry  It  is.  however,  likely  that  the  scalar  >p,^e 
required  will  be  of  order  2d  species  even  for  comparatively  simple  fuels  te.g  Lmdstedl  el  a!  [31]  and  Pope  a  a! 
[32])  and.  accordingly,  the  computational  requirements  are  exceptional.  The  current  finding  is  interesting  as  it 
suggests  that  improved  closures  for  turbulence-chemistry  interactions  combined  with  the  relative  economy  of 
unsteady  RAN*  calculations  may  provide  .  rente  towards  the  modelling  of  obstacle  wall  enhanced  turbulence 
induced  transition  to  detonation. 
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APPENDIX  \ 

The  formulations  resulting  from  she  (it. VI  [19.20]  for  the  Reynolds  stresses  and  scalar  fluxes  are  gnen  above. 
I  he  G  tensor  It-  the  following  form- 

^  A.  +  ct  jb  y  h\  — 
k  to; 

where 

=  ft#  A  *  ft*V>  +  Px 6, ,8 
i,+rAft  ~?Ah..i 
yA A,  +r<#A  *r«A  A 


The  model  contains  ’  1  coefficients  of  which  6  can  be  eliminated  by  exact  constraints  deduced  form  the  Navicr- 
Stokes  equations  The  remaining  5  were  assigned  values  hv  matching  experimental  data  f  1 9]  and  no  adjustment 
made  in  the  present  work. 
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Figure  1.  Comparison  of  computed  and  measured  flame  images  around  the  '.•.me  the  flame  moves  past  the 
obstacle.  The  computation  shows  a  later  time  and  the  schlicren  images  illustrate  the  process  leading  up  to  the 
formation  of  the  large  scale  structures 
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figure  2  Comparison  ot  experimental  anti  computed  flame  images  a!  Imer  times  3oth  the  computations  and  the 
measurements  dearly  show  large  scale  instabilities  propagating  downstream  of  the  obsiade.  The  experimental 
images  .shots  the  superimposed  turbulent  flame  wrinkling  The  large  scale  instabilities  appear  to  be  of  order  10  ’ 
m  and  the  flame  wrinkling  of  order  10  '  m 
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Station  4  (0.03  m) 


Station  2  (0.50  m) 


Figure  3  Comparison  of  computed  and  measured  pressure  traces  For  a  turbulent  gaseous  explosion  in  a  pre¬ 
existing  turbulence  field.  The  figures  show  the  pressure  traces  on  both  sides  of  the  obstacle. 
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Station  1  (0.415  m) 


Figure  -1  Comparison  of  computed  and  measured  axial  mean  velocities  at  4  measuring  stations  aiong  the  415 
mm  plane  (directh  above  the  obstacle)  for  a  turbulent  gaseous  explosion  in  a  pre-existing  turbulence  field 

Station  1  (0.415  m) 


Figure  5.  Comparison  of  computed  and  measured  cross  stream  mean  velocities  at  -  measuring  stations  along  the 
415  nun  plane  (directly  above  the  obstacle)  for  a  turbulent  gaseous  explosion  in  a  pre-existing  turbulence  field 


Station  3  (0.51 5  m) 


Time  [s] 

Figure  6.  Comparison  of  computed  and  measured  axial  mean  velocities  aT  4  measuring  stations  along  the  515 
mm  plane  i  100  mm  downstream  of  the  obstacle)  for  a  turbulent  gaseous  e\n  sion  in  a  pre-existing  turbulence 
field. 

Station  3  (0.515  m) 


Figure  7.  Comparison  of  computed  and  measured  cross  stream  mean  velocities  a:  4  measuring  stations  along  the 
515  mm  plane  (1<k)  mm  downstream  of  the  obstacle)  for  a  turbulent  gaseous  explosion  in  a  pre-existing 
turbulence  field 
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Station  3  (0.515  m) 


Figure  J>  Companion  or  computed  and  measured  axial  turbulence  veli-cities  at  4  measuring  station?  alone  the 
515  mm  plane  (  100  mm  downstream  of  the  obstacle)  for  a  turbulent  gaseous  explosion  in  a  pre-existing 
turbulence  field. 

Station  3  (0.515  m) 


Time  [s] 

Figure  9.  Comparison  of  computed  and  measured  cross  stream  turbulence  velocities  at  4  measuring  station-, 
along  the  515  mm  plane  <100  mm  downstream  of  the  obstacle)  for  a  turbulent  gaseous  explosion  in  a  pre¬ 
existing  turbulence  field. 
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